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It is widely believed that the assumption of homogeneity is a good zeroth order approximation 
for the expansion of our Universe. We analyze the correction due to subhorizon inhomogeneous 
gravitational fields. While at early times this contribution (which may act as a negative pressure 
component) is perturbatively subdominant, we show that the perturbative series is likely to diverge 
at redshift of order 1, due to the growth of perturbations. So, the homogeneous Friedmann equation 
can not be trusted at late times. We suggest that the puzzling observations of a present acceleration 
of the Universe, may just be due to the unjustihed use of the Friedmann equation and not to the 
presence of a Dark Energy component. This would completely solve the coincidence problem. 


PACS numbers: 98.80.cq 

INTRODUCTION 

Since the birth of modern cosmology the theoretical 
framework used to describe the evolution of the observed 
Universe has been based on the Einstein equations ap¬ 
plied to a perfectly homogeneous and isotropic model 
filled with various types of fluid components. It is uni¬ 
versally believed, since the Universe on large scales looks 
almost homogeneous, that the assumption of homogene¬ 
ity is a good zeroth order approximation for describing 
our Universe, and on top of this one usually adds lin¬ 
ear fluctuations (with initial amplitude of order 10“®, 
as measured by CMB experiments) in order to describe 
the growth of density perturbations and the formation 
of structures. Some steps beyond this paradigm have 
been made for describing late time behaviour of small 
scale density perturbations (which at late times become 
non-linear) and more recently in the second order rel¬ 
ativistic perturbation theory to catch new effects (as 
non-gaussianities in the spectrum of primordial pertur¬ 
bations). 

However, few efforts have been made to understand 
how reliable is the assumption of homogeneity. Of course 
there are corrections to this description due to pertur¬ 
bations: the nonlinearity of Einstein equations implies 
that the presence of perturbations modifies the back¬ 
ground. The usual way of thinking implies that these 
corrections should be small: since they are nonlinear, 
the effect should appear as a second order effect, so they 
should be of order 10“^°. 

However this reasoning is too naive, because it does 
not take into account the fact that perturbations grow 
during the matter era. 

Also, in recent years the exploration of the universe at 
redshifts of order unity has provided information about 
the time evolution of the expansion rate of the universe. 
Observations indicate a very unexpected fact: the uni¬ 
verse is mesently undergoing a phase of accelerated ex¬ 
pansion [jj]. The accelerated expansion is usually inter¬ 
preted as evidence for a negative pressure component to 


the mass-energy density of this homogeneous model uni¬ 
verse. The nature of this “Dark Energy” fluid is a huge 
mystery. 

The goal of this paper is to show that this interpre¬ 
tation can not be trusted, since the approximation of 
perfect homogeneous background becomes inadequate at 
late times due to the growth of perturbations and it fails 
at sufficiently late times. In other words the average de¬ 
scription of our Universe using Friedmann equation at 
small redshift is likely to be incorrect, due to the effect 
of the inhomogeneities. 

It is striking the fact that the Friedmann equation fails 
around the same epoch at which the acceleration is ob¬ 
served. So a Universe filled only with matter, may not 
necessarily be in contrast with a late time acceleration: 
it might just be due to the unjustified use of the Fried¬ 
mann equation. This would completely answer to the 
coincidence problem: at some time around structure for¬ 
mation the fluctuations become big enough to influence 
the background, even irrespectively of the initial ampli¬ 
tude of the fluctuations ^. 

The fact that we are pointing out here is entirely due to 
subhorizon perturbations, and so it is completely differ¬ 
ent from what has been suggested in . In that work it 
has been suggested that a large variance due to superhori¬ 
zon modes can have an important effect on the observed 
expansion: if we live in a large void we might experience 
an expansion very different from the usual FRW. 

We analyze here instead the mean effect of inho- 
mogeities on the emansion rate at perturbative level (us¬ 
ing the results of [3) and estimate the mean correction 
which is unavoidably present in any region of the Uni¬ 
verse. This may be interpreted as the contribution of 
inhomogeneous gravitational fields to the energy content 
of the Universe. It is in fact a pure nonlinear effect, due 


^ We thank N. Shuhmaher for pointing out this fact. 
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to the fact that gravitational fields gravitate.^. We will 
analyze the contribution of the subhorizon modes, ignor¬ 
ing the influence of the unknown superhorizon modes. 
We stress so that there are no free parameters here and 
in principle everything is calculable and predictable. 

To be able to catch this nonlinear effect, the first non¬ 
trivial order is second order. At this level, we found that 
the mean contribution of the second order is 0(10“®)^. 
It can also be shown that this contribution may act as an 
effective negative pressure component. This gives rise to 
a first crucial observation: the effect is five orders of mag¬ 
nitude bigger than what is naively expected, and this is 
due to the growth of perturbations after matter-radiation 
equality. 

We point out here that although this contribution may 
still appear small, one has to consider the full perturba¬ 
tive series. By a simple order of magnitude estimate we 
show that in most of the history of the Universe, the 
magnitude of higher orders is decreasing, and so the per¬ 
turbative approach is likely to give a good description. 
However, the main result of this paper is that the expan¬ 
sion parameter becomes of order 1 at redshifts of order 
z ^ 1, and not 10“® as one might naively think. This 
tells us that all the terms in the perturbative series are 
going to be of the same order at late times. There is some 
uncertainty on what is exactly the redshift at which the 
expansion parameter becomes 1, but it is certainly close 
to present time. As we said, at the very same redshifts 
the Universe is observed to accelerate. We speculate that 
the acceleration might just be the result of the fact that 
the series diverges, and that the present ACDM concor¬ 
dance model might just be regarded as an effective theory 
in which all the terms of the series have been resummed, 
giving rise to an effective negative pressure component. 


THE MEAN EFFECT OF INHOMOGENEITIES 


Newton constant and H is the Hubble parameter), and 
compute perturbatively the correction. However we are 
speculating that the correction itself drives the critical 
density pc away from pM due to the fact that H does 
not evolve as in a matter dominated homogeneous uni¬ 
verse, at late times. So, when using the expression of the 
matter power spectrum (shape factor and normalization) 
for estimating the effect, we will not use as a present 
value U® = 1'*, but a smaller value, as measured by 
observations on the turnover point of the matter power 
spectrum. On the other hand, when solving Einstein 
equations and getting the time evolution of the growth 
function for density perturbations we assume, in a per¬ 
turbatively consistent way, matter domination, and this 
will be valid until the breakdown of the approximation 
(that is z of order 1). 

We will call r the (conformal) time in this gauge, and 

the spatial coordinates, so that the metric has the form 

ds^ = a^ij) [—-|- x^)dx^dx^] ( 1 ) 

A prime will denote derivatives with respect to conformal 
time. We also introduce a proper time by dt = adr. The 
metric will be expanded perturbatively ([^) expressing 
all the metric variables as a function of a single initial 
condition , given by the function ip{x), which is constant 
in time and related to the first-order density pertur¬ 
bation, by the Poisson equation, 

vV = (2) 

We assume initial conditions as given by inflation, ne¬ 
glecting initial vector and tensor modes (|2|)- 

Then we will be interested in the local expansion rate 
0 = as seen by a comoving observer. We are going 

to average this quantity over a finite region of size R at 
fixed time r with the following prescription: 


We base our analysis on the results of 0 , and we recall 
the notations used there. We consider a Universe filled 
only by irrotational dust and choose the coordinates of 
an observer at rest with the dust (be., comoving coor¬ 
dinates), and with the same time coordinate for every 
point of the spatial hypersurface (be., synchronous coor¬ 
dinates). This system of coordinates can be chosen if the 
Universe is filled by a single pressureless component. 

Our approach is to start with a flat matter dominated 
Universe (so with an initial matter energy density frac¬ 
tion Um = ^ — = 1; where pM is the matter 

energy density, pc is the critical energy density, G is the 


V'liG x^) x'-) 

= -——- - 

v'7(t, x^) 


(3) 


After having computed (0), we compare this with the 

naive expectation = 3ib = {9) — {59)) that 

an observer would have measuring the average energy 
density {p) and applying the homogeneous Friedmann 
equations. So what we want to estimate is the quantity: 


(W) ^ (9) - 3 ^ 

3i7 3i7 

Also, we stress that here we are interested just in the 
mean corrections, which means that these are effects that 


^ This effects are often called “backreaction” in previous literature. 

^ The same result was found in [3, as an order of magnitude esti¬ 
mate. ^ With the superscript we always mean present time. 
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are present in any realization of the region of radius R. So 
this means that the extension of the r^ion R is irrelevant 


for this discussion. As we found in the mean correc¬ 
tion is due to modes which are subhorizon, and in fact the 
main contribution comes from wavelenghts of 0(20 — 50) 
Mpc, and it is insensitive to the infrared and to the ul¬ 
traviolet. These variances around this mean values are 
also interestin^as they seem to depend on superhorizon 
modes (see p,H), but we do not discuss them here and 
simply assume they are negligible (this is certainly true 
for example if fluctuations with wavelenght bigger than 
the present Hubble radius are negligible). 

The result that we found consistently at second order 


m 1^ 

1 was: 

(SO) 

1 

3H 



r2 
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-5(vV>, + 5(<*>vV> + (2{^,^0) + 


^(((vV)’)-<»■•*>,«>) 

(vV>,-|J<vV),(vV>, 


(4) 


where (...)i = f (...)d^x/(f cfix). Let us comment on 
the various terms that are present here. First, note the 
absence of terms that go as (p'^, without gradients. ® 

Then, as long as we are interested only in the statis¬ 
tical mean, many terms are irrelevant. Indeed the first 
order term has obviously zero mean, and the terms in 
the last line are small if we perform the spatial averages 
on sufficiently large scales (for example at horizon scales, 
they are completely irrelevant). 

So, only two types of terms are potentially impor¬ 
tant: the terms with two spatial derivatives + 

yI (^ip, and the terms with four spatial derivatives. 

As already noted in 113 , a term with four derivatives 
would give a correction of order 1 to the expansion at 
present times, already at second order. In fact, as we are 
going to show again, the presence of gradients boosts the 
corrections: the more gradients we have, the bigger is the 
correction. However, an important property to keep in 
mind is that the mean of a total spatial gradient is zero 
(this can be easily verified going to Fourier space (fol¬ 
lowing i) ) In other words, as long as we are interested 
only in the mean, we are allowed to derive by parts our 
expression. 

So, doing that, we immediately see that the four 
derivative terms gives zero, and that the result can be 
written as: 


m 

3H 


25 I 

Mo2iJ2 




(5) 


Now let us comment on the two derivative term. Sym¬ 
bolically we will call it term. This contribution was 
found to be of order I0“®. We stress now that this is not 
a small correction, in the sense that naively one could ex¬ 
pect a second order term to be of order I0“^°. This is the 
crucial point: the presence of the two powers of k boosts 
the contribution of five orders of magnitude approxima- 
tively at the present epoch. The reason for this is that 
the scales fco which have had more time to grow are scales 
much smaller than the horizon, and so a factor (ko/H)^ 
gives a big boost. We are going to show in detail why 
this happens, and show that this signals the breakdown 
of the perturbative expansion. In fact, extending the es¬ 
timates to analoguous terms that will appear at higher 
order in perturbation theory, we find that the expansion 
parameter is of order I at late times. 

Let us have a look at higher order terms. Of course, 
being interested only in the mean, and assuming the vari¬ 
able p gaussianly distributed, only the even orders of 
perturbation are relevant. So the next relevant order is 
the fourth order. For each order the maximum number 
of spatial derivatives grows by two. As we said the more 
derivatives we have the bigger is the contribution to the 
mean. So at fourth order the highest derivative term is 
k^p"^. As we will show, it is possible to see that such a 
term again yields a order I possible correction. At this 
point one might be tempted to speculate that at fourth 
order this may not be zero (as in 0). However, the 
fact that the highest derivative term is a total spatial 
gradient at second order may not be accidental, and one 
might expect that at any order in perturbation theory 
the term with the highest number of derivatives always 
appear as a spatial gradient. So, in order to be conser¬ 
vative, we will not rely on this speculation, and assume 
that all the terms with the highest number of derivatives 
are zero at any order. Also, as we discuss in sectionp the 
highest derivative terms are the newtonian terms and it 
is expected to get zero mean corrections at the newtonian 
level. 

Now, the next term at fourth order is k^p^, and by 
analogy with the second order term k'^p'^, we expect it 
to be nonzero. So, generically at order n, we assume 
terms like k'^’^p^ to be zero, while we keep terms like 

pn-2pn_ 

For comparison with , we stress again that here we 
are not relying on effects due to UV divergencies (that is 
to non linear structures), since our considerations involve 
effects due to scales which are still linear today. And we 
do not rely on newtonian terms, but our effect is purely 
relativistic (see section,^. 


THE EXPANSION PARAMETER 


^ This was found first in , where the expansion (without gradi¬ 
ents) was computed in the longitudinal gauge, going then back 
to the comoving coordinates. 


So let us analyze the magnitude of the k'^p'^ term. As¬ 
suming the usual flat spectrum for p, its value is given 
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by (i): 


—9 *2 ^ -/?r 

fc2(^2 ^ - 

ao 


\hMpc~^ 


dqqT^q) (6) 


where A is the amplitude of the spectrum of ip and it 
is given by 1.9 x where we took into account 

of the present value of Hm in the normalization l|l2|'l. 
Then, T{q) is the Bardeen, Bond, Kaiser, Szalay (BBKS) 
transfer function 


Tiq) = X 


[ 1 - 


2.34g 

3.89g + (16.1g)2 + {hMqf + (6.71g)^] (7) 


which describes the stagnation effect of the small scale 
perturbations that entered the horizon during the radia¬ 
tion era. Of course at small q, T‘^{q) —> 1, while at large 
q, T‘^{q) q~'^\T?q. The integral receives the biggest 

contribution from the scale at which the integrand has a 
maximum, which corresponds to 0(20 — 50) Mpc. Here 
r is the shape factor, defined for a fiat universe in terms 
of the baryon fraction Q.b and the total value of 
as r = 0^^/iexp(—Os — ). T is sensitive 

especially to the value of and it is decreasing as 

increases. Finally, the integral in eq.® has the 
numerical value of 2.3 x 10“^. 

For comparison we may estimate a term (al¬ 

though as we said it is not present in the result given 
by eq.0): 


larger contribution than the integral that we found at 
second order, and might be sensitive to the small scales. 
However, in order to make a conservative estimate we 
do not rely on such speculations and we just assume the 
minimal value 

pOO 

Int ~ / dqqT^{q) K, 2.i X . 

JQ 

At the generic order n (where n is even) so we have: 

where e is the expansion parameter: 

A f HoT Y ^ ^ 

e = ■;- TT7 - T X 

1 -I- 2 \ hMpc 1 ) 

The crucial fact is that this expansion parameter becomes 
of order 1 for small z (the value depends on uncertainties 
on the cosmological parameters, mainly through F): 


while already for z > 1, the terms of the series start 
decreasing, so one may assume that the series converges. 
The expansion parameter is plotted in fig^ for different 
values of the cosmological parameters. 

In a similar way one can also see that the highest k 
terms (that we set to zero) are of the same order of k'^ip'^ 
at small redshifts. 


k^ip^ = A^ 



HT 

hMpc~^ 


dqq^T^q) ( 8 ) 


INHOMOGENEITIES AND ACCELERATION 


Such a result would imply a big correction already at 
second order. In fact the factor large 

as ten or eleven orders of magnitude (depending on the 
value of cosmological parameters). In addition, the in¬ 
tegral is even logaritmically divergent in the ultraviolet, 
and so gets its main contribution from the small scales. 
Even putting a cutoff at the nonlinear scales, the inte¬ 
gral results in a 0(1) factor. And finally this correction 
increases with time twice as fast as the k'^p^ term, be¬ 
coming big at late times. 

Now let us go to estimate the higher orders. As we 
said, going to fourth order, we safely assume that the 
leading term is k^pd. This is estimated as follows: 





HT 


i^'hMpc~ 


6 

X Int^ 


Here the term Int^ stands for one or more integrals in 
dq. Depending on what is the precise structure of the 
contraction of indices in the k^p'^ term, it may contain 
also ultraviolet divergent integrals which would give a 


Now, we may speculate that the fact that the expan¬ 
sion parameter is of order 1 just around 2 « 1 is related 
to the appearance of acceleration in the observed Uni¬ 
verse. If this is the case, this fully solves the coincidence 
problem. Even more: one can see that the time at which 
e is of order 1 depends on the amplitude of the initial 
fluctuations, so it may appear a coincidence that the ini¬ 
tial value of 10“® is such that the expansion breaks down 
right today. However for any initial amplitude the time 
at which the series breaks down is always around the time 
of structure formation, and since we can only live around 
that time, there is no coincidence at all. 

At this point, so, let us assume that the Universe starts 
accelerating as a result of the divergence of the series. 
This, in turn, stops the growth of the density perturba¬ 
tions, and so the expansion parameter itself is prevented 
to grow and it stabilizes around the value that it takes 
when the effect starts to be important (as one can see in 

fig- ©)• 

It has been pointed out in and [ig that a decel¬ 
eration parameter defined in terms of a local expansion 
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rate can never become negative in absence of vorticity 
and in absence of negative pressure fluids. However this 
result does not apply here, since the expansion we are 
considering is not 0, but is {9). When trying to derive 
the same result to this quantity, one can see that new 
negative terms appear due to the time derivative of the 
y/j, which is present in the definition of the spatial av¬ 
erage. 

In order to show this explicitly we need to introduce 
the expansion tensor, which reads 

This can be also decomposed into its trace (which is 0 = 
and a traceless tensor (the shear tensor): 

( 11 ) 

Once we have defined the expansion rate as (0)/3, its 
deceleration parameter is given by: 

( 12 ) 

The Raychaudhuri equation gives: 

^ = -y-cr2-47rG(p + 3p) (13) 

where And where we assumed zero pressure 

and zero vorticity. Another useful relation is obtained 
taking the determinant of eqs. m and cnj, which gives 

2j dt 3 ^ ^ ■ 

Using these relations it is straightforward to obtain the 
following: 



FIG. 1: The expansion parameter e as a function of the 
redshift, for different values of the parameters (h = 0.72 ± 
0.05, = 0.14 ± 0.02, Q.bh'^ = 0.024 ± 0.001, A/Qm = 

1.9 ±0.2, as taken from WMAP onlu fll'i ])). The grey dashed 
line corresponds to the central value, while the red solid lines 
to 2<j ranges. Here we used the growth factor as in matter 
domination. For comparison the green dotted line corresponds 
to Q.M = 1. The perturbative series breaks down at small z 
when e becomes of order 1 (we also extended the plot up to neg¬ 
ative z in order to show that at some point the series badly 
diverges unavoidably.) 
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12 ^^^^ ± 


2 (v { 0 )^ ) 2 


\ w ) 


, (14) 


which is an exact equation, valid in gen eral. The first line 
corresponds to what was derived in |l5| , while the second 
line arises due to the presence of in the averages. 

This clearly shows that there are negative contribu¬ 
tions and so when the corrections overcome the back¬ 
ground value, the deceleration parameter can change 
sign. 


NEWTONIAN AND POST-NEWTONIAN TERMS 

It is very interesting to note that (see 0 ) the highest 
k terms (that we neglected) are also the terms with the 
highest power of a oc (we are using here the simple 
growth factor and background evolution of a pure matter 


FIG. 2: The expansion parameter e as a function of the 
redshift, for different values of the parameters (h = 0.72 ± 
0.05, = 0.14 ± 0.02, = 0.024 ± 0.001, A/Um = 

1.9 ±0.2, as taken from WMAP onlu fXl'j ]}). The grey dashed 
line corresponds to the central value, while the red solid lines 
to 2a ranges. Here we used the growth factor for density per¬ 
turbations of a AC DM model XlAl . in order to show that e is 
prevented to grow and it stabilizes at late times, if the Uni¬ 
verse starts accelerating. 


dominated Universe). Now, in an expansion in powers of 
the speed of light c, each power of t is accompanied by 
a power of c. So the highest k terms are the newtonian 
terms (0). 

We stress instead that the terms that we considered 
here, instead, are the post-newtonian ones. 

The following example is quite illuminating in order to 
understand why the highest k terms should give zero ef¬ 
fect. Consider a homogeneous universe filled with many 
spherical inhomogeneous regions: outside these regions 
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matter is uniformly distributed with energy density pq, 
and inside these regions matter density is inhomoge- 
neously distributed (but still maintaining the spherical 
symmetry and in such a way that the total amount of 
matter inside a sphere is equal to the amount of a same¬ 
sized sphere with homogeneous matter density po)®- 

Now, in newtonian gravity, we can use the superposi¬ 
tion principle and consider separately the effect of a single 
sphere. By Gauss’ theorem the only thing that matters 
is the total amount of matter inside the sphere, so the av¬ 
erage expansion rate over a region that contains spheres 
is the same one of a perfectly homogenous universe with 
matter density po- 

The only way in order to guarantee that the effect is 
zero in this example is that the Newtonian terms be total 
spatial gradients. In fact in this case one can rewrite the 
spatial average as a flux of a vector field over a surface. 
Taking the surface in such a way that it does not intersect 
any sphere, this gives zero because the vector field in the 
homogeneous region is zero. 

This shows that Newtonian terms have to be surface 
terms, and in order to have an effect one has to consider 
the purely relativistic terms (and the post-newtonian ap¬ 
proximation should be accurate enough). 

So, what we are claiming is that the series (in an ex¬ 
pansion in the smallness of the metric elements) of all the 
post-newtonian terms is diverging, and we may speculate 
that the resummation of all these terms may give rise to 
an accelerating Universe. Or in other words it may give 
rise to an effective theory in which the Friedmann equa¬ 
tions are still valid, but with an effective A term. In order 
to prove this kind of “phase transition” at redshift of or¬ 
der 1, so, one needs a treatment of the postnewtonian 
terms, which does not assume that the metric elements 
are small. 


CONCLUSIONS 

We have shown that, even being very conservative, the 
contribution of the inhomogeneities to the average ex¬ 
pansion rate becomes perturbatively out of control at 
low redshifts. This leads to the conclusion that the usual 
Friedmann equation cannot be trusted at late time. 

So, instead of using the naive homogeneous model to fit 
the data, one should to take a nonperturbative approach 
to Einstein equations in a inhomogeneous Universe, be¬ 
fore talking of Dark Energy. 

We suggest that the failure of this equation might be 
the reason why we need to introduce mysterious com¬ 
ponents in order to fit the data. We speculate that the 
resummation of the terms of the diverging series (which 


We thank N.Afshordi for pointing out this example. 


may act as a negative pressure component) might lead 
to an effective Friedmann equation in which the series is 
resummed, giving an effective “Dark Energy”. Once this 
effect is taken into account, the perturbations themselves 
stop growing (as happens in a AC DM model), and we 
speculate that the effect itself may be stabilized, reach¬ 
ing a fixed point, in which it might even be possible to 
predict the effective Dm ~ 0.3. 
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